#!/usr/bin/gnuplot
set encoding iso_8859_1
set term pdf enhanced color solid size 4in, 4in

#============================
#F(Q_LC), F(Q_NC), F(Q_LN)
#============================
#set xlabel "Q_{b}"
##set x2label "Free Energy (K_bT)"
#set x2tics
#set ylabel "Free Energy"
#set output "F1D.pdf"
#p 'Qb41.dat.wam2' u 1:2 w l ti "Q_{Aspe}"# , 'QPLU.smth.wam2' u 2:1 w l axis x2y1 lw 4 lt 7 ti "F({Q_{Aspe}})"
#============================
#F(Q_LC,Q_NC), F(Q_LC,Q_ATP), F(Q_NC,Q_AMP)
#============================
#set pm3d at bs
#set palette color
#set palette model CMY
#set palette defined ( 0 "black", 1 "blue", 2 "cyan", 3 "green", 4 "yellow", 5 "orange", 6 "coral", 7 "red")
#set palette defined ( 0 "black", 1 "red", 2 "orange", 3 "yellow", 4 "green", 6 "blue")
#set palette defined ( 0 "black", 1 "blue", 2 "cyan", 3 "green", 4 "yellow", 5 "orange", 6 "coral", 7 "red", 8 "white")
#set palette defined ( 0 "black", 1 "blue", 2 "cyan", 3 "green", 4 "yellow", 7 "red")
#set view map
#set cblabel 'Energy Bar' textcolor lt 3
set output "RMSDtotal.pdf"
#set multiplot layout 2,2

set title "The conformational distribution of DCL1-A"
set xlabel "RMSD (nm)"
set ylabel "Probability"
set xtics 1
plot "PD005.dat" u 1:2 w l lt 3 lc rgb "red" ti "C_{Salt}=50mM",\
     "PD01.dat" w l lt 4 lc rgb "green" ti "C_{Salt}=100mM",\
     "PD015.dat" w l lt 5 lc rgb "blue" ti "C_{Salt}=150mM",\
     "PD02.dat" w l lt 6 lc rgb "purple" ti "C_{Salt}=200mM",\
     "PD03.dat" w l lt 7 lc rgb "dark-yellow" ti "C_{Salt}=300mM"

#set title "F(COM,Q_{f})"
#set xlabel "COM"
#set ylabel "Q_{f}"
#set xtics 20
#set xrange [-0.04:101.0]
#set yrange [-0.04:101.0]
##set output "F2D_COM_Qf_org.pdf"
#splot "total.dat.wam23.org" using 2:1:3 w p pt 5 ps 0.3 palette noti
#
#set title "F(Q_{f},Q_{b})"
#set xlabel "Q_{f}"
#set ylabel "Q_{b}"
#set xtics 0.2
#set xrange [-0.04:1.0]
#set yrange [-0.04:1.0]
##set output "F2D_Qf_Qb.pdf"
#splot "total.dat.wam12" using 2:1:3 w p pt 5 ps 0.3 palette noti
